Order parameter symmetry and mode coupling effects at dirty superconducting 

quantum phase transitions 
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We derive an order-parameter field theory for a quantum phase transition between a disordered 
metal and an exotic (non-s-wave) superconductor. Mode coupling effects between the order parame- 
ter and other fermionic soft modes lead to an effective long-range interaction between the anomalous 
density fluctuations which is reflected in singularities in the free energy functional. However, this 
long-range interaction is not strong enough to suppress disorder fluctuations. The asymptotic crit- 
ical region is characterized by run-away flow to large disorder. For weak coupling, this asymptotic 
region is very narrow. It is preempted by a wide crossover regime with mean-field critical behavior 
and, in the p-wave case, logarithmic corrections to scaling in all dimensions. 



I. INTRODUCTION 

Quantum phase transitions are one of the most intrigu- 
ing problems in today's condensed matter physicsiS* 3 . 
In addition to being of fundamental interest, they 
are believed to underlie a number of interesting low- 
temperature phenomena, in particular various forms of 
exotic superconductivity^*^ 

In a seminal paper ^ Hertz introduced a general scheme 
for the analysis of quantum phase transitions in itin- 
erant electronic systems. This scheme is based on the 
Landau-Ginzburg- Wilson (LGW) approach of integrat- 
ing out the fermionic degrees of freedom and deriving a 
free energy functional in terms of the order parameter 
fluctuations only. However, in recent years, it has be- 
come clear that for many transitions, there are problems 
with Hertz' scheme because besides the order parame- 
ter fluctuations, which are soft (gapless) at the critical 
point, there are additional fermionic soft modes in the 
system. These additional soft modes exist not only at 
the critical point but also in the bulk phases. They are 
related to conservation laws and/or broken symmetries 
and constitute examples of generic scale invariancei If 
the coupling between the order parameter and these ad- 
ditional soft modes is sufficiently strong it generates an 
effective long-range interaction between the order param- 
eter fluctuations. This is reflected in a non-analytic wave- 
number dependence of the vertices of the LGW theory of 
the corresponding transition&2ii£ Generically, such non- 
localities in the LGW theory will lead to non-mean-field 
critical behavior of the quantum phase transition. 

The precise influence of the mode-coupling effects on 
a quantum phase transition depends on the structure of 
the additional soft modes and their coupling to the order 
parameter. For the clean ferromagnetic transition, the 
mode-coupling effects can either lead to a fluctuation- 
induced first order transition or to non-mean-field crit- 
ical behavior pii*i^ For dirty electrons, the transitions 
is generically of second order but with highly unusual 
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Even stronger effects were found for the 



transition between a dirty metal and a conventional (s- 
wave) superconductor. 15 Here, the mode-coupling effects 
lead to a critical point with exponential scaling, i.e., the 
correlation length behaves as £ ~ e 1 /'*', where t is dis- 
tance from the quantum critical point. Based on general 
symmetries of itinerant electronic systems, it was recently 
showni^ that homogeneous (q = 0)il order parameters 
in the particle-particle (Cooper) and spin-triplet particle- 
hole channels are generically affected by mode coupling 
effects while order parameters in the particle-hole spin- 
singlet channel do allow for a local LGW theory. 

All of the above examples are quantum phase tran- 
sitions with zero angular momentum order parameters. 
The effect of mode-coupling on order-parameters with 
finite angular momentum are much less understood. 
Herbut 18 studied the d-wave superconducting quantum 
phase transition in two dimensions within Hertz' scheme 
(which is equivalent to Gorkov's mean-field theory). He 
found that the typical Cooper channel (BCS) logarithmic 
singularities are demoted to irrelevant terms by the d- 
wave symmetry. This raises the important general ques- 
tion: How does a finite order parameter angular momen- 
tum influence the coupling between order parameter and 
additional fermionic soft modes? 

In this Paper, we study this question for quantum 
phase transitions between a metal and an exotic (non- 
s-wave) superconductor in the presence of non-magnetic 
quenched disorder. These transitions are of experimen- 
tal importance since various superconducting states with 
p,d and maybe higher symmetries have been observed 
recently, and their quantum phase transitions are exper- 
imentally accessible 4,5,6 in principle. Experiments per- 
formed on the weakly ferromagnetic compounds \JGe^ 
and ZrZn^ 6 " revealed the existence of a superconducting 
phase within the ferromagnetic phase at temperatures 
below IK. It is believed 5 - that both superconductivity 
and ferromagnetism arise from the same electrons. One 
possible mechanism for this type of superconductivity is 
p-wave triplet pairing mediated by magnetic fluctuations 
in the vicinity of a magnetic quantum critical pointf^iS 3 . 
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although this is still not a settled issue. The onset of the 
phenomenon has proven to be very sensitive to the pres- 
ence of non-magnetic disorder, making it observable only 
in highly pure samples. This fact also points to a non- 
s-wave order parameter. In ZrZn2 the superconducting 
quantum phase transition as a function of disorder has 
already been observed^ 

Our results can be summarized as follows: Mode- 
coupling induced singularities exist for all order param- 
eter angular momenta L. However, with increasing L, 
they are more and more suppressed (by a factor |q| 2i .) 
As a result, the LGW theory of our superconducting 
transition is equivalent to that of the itinerant antiferro- 
magnetic transitioniii^iSLSi While the asymptotic crit- 
ical behavior of this theory is not understood because 
of runaway flow to large disorder, we also show that for 
weak bare disorder the asymptotic region is exponentially 
narrow. It is preempted by a wide crossover regime with 
mean-field critical behavior and (in the p-wave case) log- 
arithmic corrections to scaling. The paper is organized as 
follows. In Scc.[nJ we derive the LGW free energy func- 
tional. In Sec. IIIII we study the LGW theory by means 
of the renormalization group and determine the critical 
behavior. In Section llVI we analyze our findings from a 
mode-coupling point of view, and discuss differences be- 
tween paramagnetic and ferromagnetic as well as clean 
and dirty systems. 



II. LANDAU-GINZBURG- WILSON THEORY 
A. p-wave pairing case 

In this section, we derive an effective LGW theory for 
the disorder-driven quantum phase transition between a 
paramagnetic metal and a p-wave triplet superconduc- 
tor. Our starting point is a microscopic action for in- 
teracting electrons in d > 2 dimensions and subject to a 
static, single-particle random potential V(x). We assume 
a Gaussian distributed potential with [V (x^V (x 2 )]dis = 
W8(xi — X2), with W being measure of disorder strength. 
The partition function Z can be written as a functional 
integral over Grassmann variables 

Z = jD[^]e s ^ . (1) 

We decompose the action S = S p + So into the p-wave 
interaction part S p and a reference system So which com- 
prises the single-particle part, the random potential and 
Sint (the interaction in all channels other than the p-wave 
channel) : 

f - V 2 

So = / (foViW[9 T + — + ^(x) + (2) 

+ I dxJ2^(x)V{x)i>Ax)+Sint, (3) 



S P = / dxh (7(7 ,(x)-n (Tl(T , i (x) . (4) 

We use a (d + l)-vector notation, with x = (x, t), k = 

Ck.fi), fdx = J v d d xJ P dr and £ fc = £ k r£ n , x is 
a real space coordinate, t imaginary time, k momentum 
vector and f2 Matsubara frequency. n CTCT / (x) is the p-wave 
anomalous density whose Fourier transform in terms of 
the fermion fields is given by 

iw(g) = E 6k tl>„(k + q/2)il> a ,(k-q/2), (5) 

k 

with e k = k/|k|, a, a' being spin indices and • denoting 
scalar product in the orbital space. Due to the Pauli prin- 
ciple the spin state of the Cooper pair has to be a triplet, 
i. e. ua' e {(H), (II), 1/V2(U + IT)}- Which combina- 
tion of the three possible triplet components is actually 
realized depends on the system under consideration. The 
reference ensemble So describes interacting electrons in 
the presence of non-magnetic quenched disorder and no 
bare interaction in the p-wave Cooper channel. (A non- 
vanishing interaction in this channel will be generated in 
perturbation theory.) So thus describes a general system 
of disordered interacting electrons with the only restric- 
tion being that it must not undergo a phase transition in 
the parameter region we are interested in. 

A standard procedure— is used to derive a LGW order- 
parameter field theory. We decouple the interaction term 
using a Hubbard-Stratonovich transformation 24-25 intro- 
ducing a complex field (x) which plays the role of 
an order parameter. Quenched disorder treated using 
the replica trickj2& and fermionic degrees are integrated 
out, leading to an expression for the critical part of the 
partition function in terms of order-parameter only: 

Z = J D[A] e~*[ A l . (6) 

Since our emphasis is on the mode-coupling effects, and 
in order to avoid unnecessary complications in notation, 
we restrict our analysis to a certain spin component 
((jo-' = (IT)) 01 the order parameter. The LGW free en- 
ergy $[A] is expanded in powers of the order parameter 
A = A||. Up to quartic order it reads: 

$[A] = ^A"( q )(i-r tX ( 2 )( (? ))A"(g) 

- E T 2 tX ^( qi , q2 , q3 )A a ( qi )A a ( q2 ) 

91--83 

x A' 3 ( (Z3 )A^(( 7l +( 7 3-92) + 0(A 6 ), (7) 

where t]^ = T t . Here are replica indexes. The co- 
efficients of the LGW functional are determined by the 
2-point and 4-point anomalous density correlation func- 
tions of the reference system So which can be written as 
%( 2 ) = (n a n a )o and Xal = (n a n /3 n t3 n a }o (with the spin 
and component indices suppressed). 
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FIG. 1: Contributions to the leading terms of the Gaussian 
part of the LGW functional, (a) provides a constant (~ Nf) 
and the frequency dependence, \uj\t, while (b) gives the lead- 
ing momentum dependence, |q| 2 log(l/|q|). 



B. Anomalous density correlation functions in the 
p-wave channel 

In this subsection we use diagrammatic perturbation 
theory in disorder and interactions to calculate x an d 
X^ 4 -* of the reference ensemble So, focusing on the be- 
havior for weak disorder. Thus, we neglect all diagrams 
with crossed impurity lines, i.e., all weak localization ef- 
fects, which become important only at higher impurity 
concentrations. This is justified for the ferromagnetic 
superconductors*^, where the superconducting quantum 
phase transition occurs at very small disorder. 

We start our analysis by examining the 2-point correla- 
tion function, 

x (2) _ (fja n <*) . 

The leading contributions 
are obtained from the diagrams shown in Fig. ^ (details of 
the calculation of these diagrams are given in Appendix 
IX}) . Here, the external vertices represent anomalous p- 
wave densities, the solid lines are fermionic propagators 
in Born approximation, 

G- 1 (k,w) = io;-e k , <T + /i+(i/2T)sgn(a;), (8) 

where e^o- is the dispersion relation and r is the scatter- 
ing time. The double line represents a particle-particle 
impurity ladder: 



W R (q,Q,u)=W 



1 i{fl(fl + uj)<0 



( 9 ) 

£ = kpr/m is the elastic mean free path and W = 
1/(2itNft) with Np being the density of states at the 
Fermi level. The calculation of the diagrams in Fig. 
for |q| = lu = leads to X (2) = (N F /3) ln(2e F r). The 
well-known logarithmic Cooper channel (BCS) singular- 
ities are cut-off by the disorder, reflecting the suppres- 
sion of exotic superconductivity by non-magnetic scat- 
terers in analogy to the suppression of s-wave supercon- 
ductivity by magnetic impurities. 27 ' 28 [We note that, in 
contrast, s-wave superconductivity is not influenced by 
weak non-magnetic scatterers, as is signified by Ander- 
son's theorem. 29 ] However, a closer investigation of dia- 
gram^i) for finite |q| reveals that a non-analyticity of the 
form |q| 2 ln(l/|q|) survives. Thus, the p-wave symmetry 
has demoted the BCS singularity to quadratic order in |q| 
because each of the renormalized external vertices picks 
up an extra power of |q|^ 




k+q 

C2+CQ 



FIG. 2: Leading singular contributions to Xdiag' After expan- 
sion in small q, the leading order terms of (a) and (b) cancel 
(see Appendix let . 



In addition to the BCS logarithms x^ contains non- 
analyticities similar to that in the itinerant ferromag- 
net. They are caused by the leading corrections to scal- 
ing at the dirty Fermi liquid fixed point^i and can be 
viewed as particle-particle analogs of the well known 
Altshuler-Aronov corrections to density of states and 
conductivity^ For s-wave order parameters these sin- 
gularities (which only arise for interacting electrons) are 
of the form |q| d ~ 2 iili For p-wave order parameters, an 
inspection of the corresponding contributions (details 
see Appendix [5J) reveals that they are suppressed by 
a factor |q| 2 by the same mechanism as the BCS log- 
arithms. [Note that an analogous suppression occurs in 
the particle-hole channel, as can be seen from a power 
counting analysis of the Altshuler-Aronov correction to 
the conductivity: 8a = 5(jj)/cu ~ cj( d ~ 2 )/ 2 ~ |q| d ~ 2 . 
Thus, the correction to the current-current correlation 
function (jj), which is proportional to the p-wave den- 
sity, behaves as 5{jj) ~ |q| d . This means it has picked 
up an additional factor |q| 2 compared to the zero angular 
momentum channel.] 

As a result, we find that the leading singularities in 
%( 2 ) are the Cooper channel logarithms, and the leading 
behavior of x^ for (|q|,£j) — * is given by: 

= ^[ln(2 er .r)-MT-^]{ (j + 

where i,j are the order parameter component indices. 
The anisotropic q-dependence in the last term reflects 
the spatial anisotropy of the p-wave order parameter. 

We now turn our attention to the 4-point correlation 
function, x^which can be split into a replica-diagonal 



part and a replica-off-diagonal part, x 

,(4) 



(-D 
a/3 



°aPXdiai 



Xoff ^ detailed discussion of our calculation is given in 
the Appendix |0 We find that the leading contributions 
to Xdiag m t ne long-wavelength, low-frequency limit are 
coming from the diagrams shown in Fig. [21 While each 
of the two diagrams individually diverges for (q, to) — > 0, 
their leading singularities cancel, and the remaining con- 
tribution is finite and proportional to NpT 2 . The lead- 
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FIG. 3: Leading contribution to the X^off ( see Appendix ICl 



ing contribution to the replica-off-diagonal part of x^ 1S 
produced by the diagram shown on the Fig. Thus, we 
finally obtain: 



(4) 
X off 



A (N 2 /k F ) F (cu,w 4 



(11) 
(12) 



where and F a are dimensionless functions with val- 
ues between and 1 and Ad and A a are dimensionless 
prcfactors of order one. 

The results 1)110 and 1)120 have been obtained from low 
order perturbation theory. Within perturbation theory, 
it is non-trivial to prove the leading non-analyticity to 
all orders. Therefore we follow the guidance of the cor- 
responding results for the s-wave case which have been 
rigorously established using Q-field-theory and renormal- 
ization group arguments. 31 (Corresponding rigorous re- 
sults for finite angular momentum modes do not yet ex- 
ist.) Indeed, simply correcting the s-wave results from 
Ref. for the g-dependence of the renormalized ver- 



tex leads to x 



2)1 



|4-2n 



in agreement with 1 )110 and 



1)120 . Note that the singularity becomes stronger in the 
higher order anomalous density correlation functions, in 
agreement with general mode-coupling arguments^^ 3 - 

Inserting eqs. l (TDjl - lfP2j l into Eq. Q we obtain the 
LGW functional up to quartic order in A, 



$ = 



^A™(< Z )(t + |c i ;|+ C2 |q| 2 )A-(( ? ) 

q,a 

J2A*«t( q )c imj [-l + iln(^)]A^(,) 

ij 

U^2 f dvdr |A Q (r,r)| 4 

J drdrdr' |A tt (r,T)| 2 |A /3 (r,r')| 2 . (13) 



a/3 



Here we have scaled the order parameter with 
(TtNpT) 1 / 2 and replaced the quartic coefficients by num- 
bers which is sufficient for power counting purposes. The 

coefficients are: t ~ 1/r 1/T — ln(2ejrr) 



C2 



Cl 



£ 2 /t, U ~ 1/Np, and V ~ 1/ (kpr 2 ) , with f = T t N F be- 
ing a dimensionless measure of the interaction strength. 
The parameter t represents the distance from the quan- 
tum critical point. Generically, U > 0, which leads to a 
second order transition. — This completes the derivation 
of the LGW theory. 



C. Higher angular momentum channels 

In this subsection we generalize the findings to pair- 
ings in higher angular momentum (L) channels. For an- 
gular momentum L > 0, the renormalized anomalous 
density vertex is proportional to |q| L . The leading non- 
analyticity in the static anomalous density susceptibility 

Xi (l) is given by the BCS logarithms in diagram \\]p. 
They take the form 



*xi 2) (q) 



%f (<?)<(<?)) 



q|^log(l/|q|), (14) 



i.e., they are suppressed by a factor |q| 2i compared to 
the s-wave case. Here njf(q) — ^ k Yl ! {e^)ip a {k + 
q/2)tJ> rT i(k — ?/2) is a component of the anomalous den- 
sity for angular momentum L. Note that for L > 1 the 
BCS logarithm is subleading compared to the regular |q| 2 
term coming from diagram \T£l) while in the p-wave case 
the BCS logarithm provides the leading wave-number de- 
pendence in the LGW functional. The same mechanism 
also suppresses the interaction induced mode-coupling 
singularities related to corrections to scaling at the dirty 
Fermi liquid fixed point. An explicit calculation outlined 
in Appendix [B] shows that these mode coupling singular- 
ities behave at most like |q| 2L |q| d ~ 2 (d is the spatial di- 
mensionality). Therefore they are sub- leading compared 
to the BCS logarithms for all d > 2. We now turn to 
the 4-point anomalous density correlation function Xl ■ 
Because of the momentum dependence of the renormal- 
ized anomalous density vertex, xt, ' picks up an extra 
power of |q| 4i compared to the s-wave case. More gener- 
ally, any 2n-point anomalous density correlation function 
^,(2n) pj c k s U p an extra power of |q| 2ni compared to the 
s-wave case, i.e. x^ 2n ^ ~ |q| 4 +( 2i ~ 4 )". Therefore their 
singular contributions are demoted to sub-leading order 
and do not play a role for the critical behavior. 

As a result, the leading terms in the LGW functional 
for d-wave and higher order parameter angular momen- 
tum take the same form (01 30) as in the p-wave case except 
for the missing logarithmic wavenumber dependence in 
the Gaussian parti 1 ^ 



III. RENORMALIZATION GROUP ANALYSIS 

In this section, we analyze the LGW theory 1)130 by 
means of the renormalization group. There is a Gaus- 
sian fixed point with mean-field static exponents v = 1/2, 
7 = 1, r\ = 0, and a dynamical exponent of z = 2. In 
the p-wave case there are logarithmic corrections to the 
mean- field behavior in all dimensions. In order to check 
the stability of this Gaussian fixed point we study the im- 
portance of quantum and disorder fluctuations. The scale 
dimensions of U and V at the Gaussian fixed point can 
be calculated by power counting. We obtain [U] = 2 — d 
and [V] — 4 — d. Thus, the conventional fluctuation 
term (the U term) is renormalization group irrelevant 
for d > 2, but the disorder term (the replica-off-diagonal 
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Disorder starts to Gaussian fixed point 
physics 




FIG. 4: Schematic flow diagram on the critical surface. The 
Gaussian fixed point (G) is unstable; the flow goes toward 
large disorder V. For weak bare disorder, the flow stays close 
to the [/—axis until it almost reaches the Gaussian fixed point 
before crossing over (black dots) to the asymptotic destina- 
tion. The dashed line separates the region described by the 
Gaussian fixed point from the strong disorder region. 



order coupling constant Vr — V/yte^ ~ 1. Now, the 
quantum phase transition occurs at t = which leads to 
l/r c = ln(2ei?r), with T c being the dimensionless critical 
coupling. Thus, the quantum phase transition occurs at 
an exponentially small bare disorder strength which im- 
plies that Vr ~ 1 requires an exponentially large length 
scale. Using the results above we find a Ginzburg type 
criterion: 



|r-r c 



< exp 



1 



(15) 



Therefore, disorder effects become important only inside 
an exponentially narrow region surrounding r c . This 
asymptotic critical region is preempted by a wide Gaus- 
sian crossover region (region below dashed dashed line 
on Fig. ^ with mean-field critical behavior. For p-wave 
symmetry there are logarithmic correction to the power- 
law scaling whose the |q| dependence reflects the under- 
lying p-wave symmetry of the order parameter. 



quartic V term) is relevant for d < 4. In three dimensions 
the Gaussian fixed point is unstable with respect to the 
disorder term, and thus the calculation of loops is neces- 
sary to determine the asymptotic critical behavior. This 
includes the possibility of replica-symmetry breaking in 
the replica-off-diagonal quartic term. 

Rather than carrying out this program explicitly, we 
use the analogy between our transition and the disor- 
dered itinerant antiferromagnetic quantum phase transi- 
tion to discuss the asymptotic critical behavior: Except 
for the logarithmic corrections, the LGW theory, eq. i|13|) 
is identical to that of a disordered itinerant antiferromag- 
nct. This transition has been investigated in great detail 
in recent years ^22*21 By taking into account rare disor- 
der fluctuations it was found that there is no critical fixed 
point in the perturbative region of parameter space, and 
the asymptotic critical behavior is characterized by run- 
away flow toward large disorder (Fig. [3} rendering the 
perturbation expansion unjustified. The physical impli- 
cations of this runaway flow are not fully understood so 
far. Possible scenarios include a non-perturbative fixed 
point with conventional (power-law) scaling, an infinite 
randomness fixed point (relative magnitude of inhomo- 
geneities increases without limit under coarse graining), 
resulting in activated scaling, or a complete destruction 
of a sharp phase transition. Thus, from the analogy with 
the quantum phase transition in itinerant antiferromag- 
nets we conclude that the asymptotic critical behavior of 
our theory is unconventional, and, at present, unknown. 

However, in many relevant experimental systems the 
bare disorder is actually very small. Thus, one may 
ask at what length scale disorder effects start to play 
a role. The crossover scale between the Gaussian and 
the asymptotic critical behavior can be determined from 
the condition that the renormalized dimensionless dis- 



IV. CONCLUSIONS 

In this paper, we have studied the quantum phase tran- 
sition from a dirty metal to an exotic superconductor. 
Starting from a microscopic action of disordered inter- 
acting electrons, we have derived the LGW theory for 
this quantum phase transition which proved to be equiv- 
alent (up to logarithmic corrections in the Gaussian part 
in the case of p-wave pairing) to the extensively stud- 
ied LGW theory of the dirty itinerant antiferromagnetic 
transition. A renormalization group analysis yielded run- 
away flow toward large disorder. As a result, the asymp- 
totic fate of the quantum phase transition is presently 
unknown. However, we could derive a Ginzburg- type 
criterion for the importance of the disorder fluctuations. 
For weak bare disorder, as is realized in many experi- 
mental systems, the true asymptotic behavior is observed 
only exponentially close to the quantum critical point. It 
is preempted by a wide region with mean-field behavior 
(and logarithmic corrections for p-wave pairing). In this 
last section we analyze our results from a general mode- 
coupling point of view, and we also discuss experiments. 

In deriving the LGW functional, we have paid particu- 
lar attention to the coupling between the order parameter 
and additional soft modes. Mode-coupling induced singu- 
larities are indeed present in all angular momentum chan- 
nels, but they are increasingly suppressed for higher an- 
gular momentum: In the static order parameter suscepti- 
bility the singular terms pick up an extra power of |q| 2L . 
This suppression can be understood as follows: In the 
presence of non-magnetic quenched disorder, the domi- 
nant electronic soft modes are those that involve fluctu- 
ations of the number density, spin density, or anomalous 
density in the zero angular momentum channel^! while 
the corresponding densities in higher angular momen- 
tum channels are not soft. Since the different angular 
momentum modes are orthogonal at zero wavenumber, 
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the coupling between a finite angular momentum order 
parameter and the zero angular momentum soft modes 
must involve powers of the wave number |q|. These ar- 
guments suggest a very general difference between the 
mode-coupling effects in clean and in dirty electronic sys- 
tems. While the only soft modes in the dirty case are in 
the zero angular momentum channel, in clean systems, 
the charge, spin, and anomalous density fluctuations in 
all angular momentum channels are soft (corresponding 
to an infinite number of Fermi liquid parameters). There- 
fore, one expects the mode coupling singularities in a 
clean system not to be suppressed by a higher order pa- 
rameter angular momentum. This is known to be true 
for the Cooper channel logarithmic singularities which do 
not pick up an extra |q| 2i in clean electronic systems. A 
systematic investigation of this question will be published 
elsewhere i££ 

The explicit calculations in this paper were for a su- 
perconducting quantum phase transition in a paramag- 
netic system. We now discuss to what extend the results 
change if the transition happens in a ferromagnetic sys- 
tem. Let us assume the magnetization being in z direc- 
tion. Obviously, not all possible order parameter compo- 
nents are equivalent. Specifically, the symmetric triplet 
l/\/2(Ti + IT) (f° r P or /-wave pairing) as well as the 
singlet l/V2(tJ. — It) (f° r s an d d-wave pairing) are 
suppressed because the exchange gap cuts off the Cooper- 
channel singularities. In contrast, for equal spin pairing 
(|| and ||), the leading behavior is the same as discussed 
in Sections ILTI and ILTTI of this paper. 

Possible candidates for an experimental verification of 
our predictions are the ferromagnetic superconductors 
UGe^ or ZrZn2i& For these systems, a likely mechanism 
for superconductivity is p-wave triplet pairing mediated 
by magnetic fluctuations due to the vicinity to a mag- 
netic quantum critical point ^ although this has not yet 
been established beyond doubt. In ZrZii2 the vanishing 
of superconductivity as a function of disorder has actu- 
ally already been observed^ A systematic study of this 
transition would therefore be very interesting. 36 
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APPENDIX A: TWO-POINT LGW VERTEX 

In this Appendix we sketch the derivation of the ex- 
pression (|10|l for the anomalous density susceptibility x^ 2 -* 



in 3d. In a suitable parametrization and for p-wave pair- 
ing, diagram Fig. can be written as 

D% (q, u)=Tj2 ^ i (e k )y i ''(e k )G(k+q, n+w)G(-k, - fi), 

( A1 ) 

with i, j being the order parameter component in- 
dices and q, uj the external momentum and frequency. 
Yff(e-k) = Yjf (9, <fi) is a spherical harmonic and Green's 
function G is given by eq. (JSJ. A straightforward calcu- 
lation leads to: 



D a 



N F 



log(2e F r) - \uj\t - 



1 

To 



^ 2 |q| 2 



l = kpT /m is elastic mean free path and Np the density 
of states at the Fermi level. Similarly, the diagram in 
Fig. IB) is: 



D 



N F 
—i 
9 



2 <li<lj l0g(wr 



£ 2 |q| 2 



(A3) 



Adding eqs. (|A2|) and (|A3|) completes derivation of eq. 
(|ll)fl . For general angular momentum L, the analogous 
calculation shows that the BCS logarithm in D b picks up 
an extra factor |q| 2L compared to the s-wave case. 



APPENDIX B: INTERACTION EFFECTS 

In this Appendix we analyze the leading corrections to 
due to the interactions Si n t in the reference ensem- 
ble 5*0. They can be understood as corrections to scaling 
at the dirty Fermi liquid fixed point 31 and are particle- 
particle analogs of the well known Altshuler-Aronov cor- 
rections to density of states and conductivity— We first 
consider a paramagnetic reference ensemble. To first or- 
der in the interactions, the relevant diagrams are those 
in Fig. and their counterparts with bare external ver- 
tices. The wiggly line represents the interaction which 
is assumed to be short ranged and can thus be approxi- 
mated by a number T aa i (where u, a 1 denote the spin at 
the two ends of the interaction line). 

Particular attention must be paid to the diagrams a, 
b and c with bare vertices. In these diagrams the spher- 
ical harmonics in the two external vertices are not inde- 
pendent. Therefore, their contributions can potentially 
produce stronger terms than |q| 2i |q| d_2 . However, it 
turns out that these contributions do not produce any 
non-analytic terms and only contribute to the regular 
terms. In the remaining contributions the angular vari- 
ables of the external vertices are independent. They can 
be analyzed along the lines presented in Ref. HI After 
straightforward calculation one finds that the interaction 
corrections produce singularities of at most of the order 
|q|2L|q|d-2 means they are suppressed by a factor 

|q| 2L compared to the zero angular momentum case^ 

The above conclusion is easily generalized to ferromag- 
netic reference ensembles: If the magnetization is in z- 
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FIG. 5: Diagrams arising in the first order perturbation the- 
ory in interaction of the reference ensemble So, and produce 
non-analytic, |q| 2L |q| d_2 , term, which is a consequence of the 
mode-coupling effects. 



direction, the |T and J, J, components of the order param- 
eter have the same type of singularity as discussed above, 
while the leading singularities in the TJ, components are 
cut-off by the exchange gap. 



APPENDIX C: 4-POINT VERTICES 

In this Appendix we present details of calculation of 
four-point susceptibility x m 3d, eqs. (Illfl and (|12|l . 

We start with replica diagonal part, Xdiag which is ca l~ 
culated from the diagrams shown on the Fig. [21 The 
most singular contributions are produced if the frequen- 
cies structure permits all four external vertices to be 
rcnormalized by an active (retarded-advanced) ladder (4- 
ladder diagrams). A direct calculation of the first dia- 
gram (Fig. [20 leads to: 



D„ 



2£ 5 r 2 |q| 4 cos 4 (a) 
8l7r 



E 



e(o(fi + w)) 



\2n + uj\T + 



i 2 M 2 



x(l - 9\2Q + u\t - £ 2 |q| 2 )(Cl) 

An analogous calculation for diagram [2)) gives (the fre- 
quency constraint requires the extra impurity line to act 
as a single impurity line rather than an active ladder): 



D h = 



^T 2 |q| 4 cos 4 (a) 



6(0(fi + cj)) 



81tt 



|2fi + uj\t + 



tn 



-) 4 



Each of these two diagrams individually has an diver- 
gence ~ l/|q| 2 - However, because the relative combina- 
tional factor of Db is 2, the divergent contributions can- 
cel, rendering the final value for Xdian fhfite (Eq. (JTTJ). 



Similar cancellations among individually diverging di- 
agrams take place in the replica off-diagonal case, with 
the strongest singularities coming from diagrams with the 
largest number of active ladders. A set of such diagrams 
is shown on Fig. Here, at most six ladders can be 
active, leading to an infrared singularity in each of the 
diagrams of the form l/|q|. Similar calculations to the 
ones carried out above reveal that the singular contribu- 
tions of diagrams (a), (b), and (c) canceled each other. 
The remaining contribution is finite and can be estimated 
from the simple diagram Fig. [3J 



We emphasize once more that all the results for the 
singularities in the anomalous density correlation func- 
tions in the Appendices A to D have been obtained in 
low-order perturbation theory. Within perturbation the- 
ory one cannot prove that the terms obtained indeed 
represent the leading singularities to all orders. We are 
nonetheless confident that we indeed identified the lead- 
ing terms, because in the s-wave case we reproduce the 
known rigorous results from Q-field theoryii^i 




x(l- 12|2S1 + c^|t — 2£ 2 |q| 2 )(.C2) 



FIG. 6: Diagrams contributing to the Xoff Ah 1 four vertices 
are renormalized, with two extra active particle-particle lad- 
ders connecting two fermionic loops. The zeroth order terms 
in small q and ui of these diagrams cancel. 
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